Thermally- and mechanically-driven quantum turbulence in helium II 
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In most experiments with superfluid helium turbulence is generated thermally (by applying a 
. . . heat flux, as in thermal counterflow) or mechanically (by stirring the liquid). By modelling the 

04 ' superfluid vortex lines as reconnecting space curves with fixed circulation, and the driving normal 

fluid as a uniform flow (for thermal counterflow) and a synthetic turbulent flow (for mechanically- 
^^ . driven turbulence), we determine the difference between thermally-driven and mechanically-driven 

Cn ' quantum turbulence. We flnd that in mechanically-driven turbulence the energy is concentrated 

{^jqI at the large scales, the spectrum obeys Kolmogorov scaling, vortex lines have large curvature, and 

i~^ . the presence of coherent vortex structures induces vortex reconnections at small angles. On the 

contrary, in thermally-driven turbulence the energy is concentrated at the mesoscales, the curvature 
is smaller, the vorticity fleld is featureless and reconnections occur at larger angles. Our results 
V^ , suggest a method to experimentally detect the presence of superfluid vortex bundles. 



I. INTRODUCTION 

Recent work [l|, y] has highlighted similarities between the turbulence of supcrfluid helium II (quantum turbulence) 
and the turbulence of ordinary (classical) fluids. In particular, experimental 0, 3 and theoretical |5l-[ll| studies have 
established that the distribution of the supcrfluid kinetic energy over the length scales (energy spectrum) obeys the 



same k ^'^ Kolniogorov scaling of ordinary turbulence 12| where k is the wavenumber. The similarity is remarkable 



because helium II is unlike an ordinary fluid: firstly, it has a two- fluid nature, consisting of a viscous normal fluid 
component and an inviscid superfluid component coupled by a mutual friction [13| ; secondly, superfluid vorticity is not 
a continuous field (like in an ordinary fluid) but is restricted to discrete vortex filaments around which the circulation 
is fixed to the ratio of Planck's constant and the mass of one helium atom. 

In most experiments turbulence in helium II is excited mechanically (by stirring the helium with grids or propellers 
0,0,13 or forcing it along pipes [lal), or thermally (by the application of a heat flux [la, [131) ■ The simplest, most 
studied form of thermal stirring (to which hereafter we restrict our work) is called thermal counterflow. The name 
arises because the normal fluid and the superfluid move in opposite directions, so that their velocity difference is 
proportional to the applied heat flux and the net mass flux is zero. Other forms of heat transfer (e.g. pure superflow) 
and other techniques to generate turbulence (e.g. ultrasound [l^, ion injection and spin-downs [l^) are either less 
studied, or refer to the low temperature limit (below 1 K), or have a special character (rotating turbulence, turbulent 
fronts, the Kibble-Zurek mechanism, etc.). thus they are not our interest here, and nor are the special methods used 
to model them; for these aspects we refer the reader to a recent review [20] . 

The aim of this work is to clarify the difference between thermally-excited counterflow turbulence and mechanically- 
excited turbulence in helium II. For simplicity, we are concerned only with statistical steady state turbulence away 
from boundaries (thus ignoring the important problems of turbulence decay and flow proflles) , and at the relatively 
high temperatures, where mutual friction plays a role. 

After setting up the necessary numerical models (section |lT| we compute for the first time the energy spectrum 
of counterflow turbulence, compare it to the spectrum of mechanically-driven turbulence, and test the idea which 
has been proposed in the literature that "counterflow turbulence has only one length scale", meaning the average 
intervortex distance fsection llll)) . We also find that thermally- and mechanically-induced turbulence differ with respect 
to curvature (section |IV|, the presence of coherent structures (section |V| and vortex reconnection statistics. The last 
result suggests a method to detect experimentally the existence of superfluid vortex bundles (section |VI|. Section IVlII 
summarises the conclusions. 



II. NUMERICAL METHOD 

Following Schwarz \2l\ . we model quantum vortex filaments as space curves s(^,t) which move according to 

iot I „f ^, /, ext tot\ I „l ^, ( „l ^, [ ext tot\\ /i \ 

— =V3 -I- as X (v„ -v^ )-as x (s x (v„ - v, )), (1) 

where t is time, a and a' are temperature dependent friction coefficients |22| . s' = ds/d£, is the unit tangent vector 
at the point s, ^ is arc length, and vjf'* and v^^* are externally applied normal fiuid and superfluid velocities. The 
self-induced velocity of the vortex filament at the point s is given by the Biot-Savart law [2j| 

v^ = -f / ii^ X dr, (2) 

47r Jc |s - r|3 

where k = 9.97 x 10"'' cm^/s is the quantum of circulation and the line integral extends over the entire vortex 
configuration C. The total superfiuid velocity is thus v*°* = v* -I- v^^* . 

In the case of thermal counterflow (which hereafter we simply refer to as thermally-driven, to contrast it to 
mechanically-driven turbulence), v^^* and v|^* arise from the imposed heat flux Q, where u„s = {vf^* — v^^*| = 
Q/{psST) is the counterflow velocity, w^^* = Q/{pST), wf^* = —pnvf{^*/ps, S is the specific entropy, T is the tempera- 
ture, ps and pn are the superfluid and normal fluid densities, and p = Ps + Pn- We make the usual [2l|,[2J] simplifying 
assumption of uniform velocity profiles v^"^* and vf^* away from boundaries, neglecting the possibility that at suffi- 
ciently large Q the normal fluid becomes turbulent [25|, [2a] . It is also convenient to perform the calculation in the 
frame of reference of the imposed superflow, setting v^^* = 0. 

In the case of mechanically-induced turbulence, since we ignore boundaries and flow proflles, we set v^^'* = in 
Eq. (ID) and replace vj^^* with the following synthetic turbulent flow [27| : 



m=M 

v^^*(s, t) = X! (^" ^ '*'" '=°'^<^'" + ^™ ^ ^™ sin0„), (3) 

771 — 1 



where 0„i = k„j -s + cjmt, k,„ and iOm — ^k'^E(km) are wave vectors and angular frequencies. This v^^* is solenoidal, 
time-dependent, and, with a suitable choice of Am and Bm (adapted to the periodic box [28|), its energy spectrum 
has Kolmogorov form E{km) ~ km in the range from ki (corresponding to the integral scale) to fcj\/ (corresponding 
to the dissipation scale). Synthetic turbulence is widely used to study transport properties, and compares very well 
with direct numerical simulations and experiments (for example it satisfies observed two-points turbulence statistics) . 

Our calculations are performed in a periodic cube of size D = O.l cm. The numerical techniques to discretize the 
vortex filaments into a variable number of points Sj (j = 1, • • ■ N) held at minimum separation A^/2, compute the 
time evolution, de-singularize the Biot-Savart integrals, evaluate v* using a tree-method (with critical opening angle 
0.4), and algorithmically perform vortex reconnections when vortex lines come sufficiently close to each other, are all 
described in our previous papers [3, [29l - [3l| . 

It must be stressed that our models have a limitation: the normal fluid is prescribed rather than computed self- 
consistently. The inclusion of the back-reaction of the superfiuid vortices onto the normal fluid would require the 
numerical solution of the Navier-Stokes equation for the normal fluid (suitably modified by the inclusion of a mutual 
friction term), alongside the time evolution of the superfiuid vortices. However, a dynamically self-consistent model 
would be very complex and computationally expensive, and one could not easily explore parameter space and the 
effects of changing numerical resolution and initial conditions. This approach was attempted for a single vortex ring 



321. In the case of turbulence, this approach has so far been limited to the initial growth of a cloud of vortex lines 
331. 



III. ENERGY SPECTRUM 

We choose temperature T = 1.9 K (at which a = 0.206 and a' = 0.0083) which is typical of experiments and allows 
direct comparison with previous work, and numerical resolution A^ = 0.0016 cm. First we calculate thermally-induced 
turbulent vortex tangles at increasing values of u„s. We find that, after an initial transient, the vortex line density 
L (vortex length per unit volume) saturates to a statistically steady state (see Fig. [1]) of density L = 7^f^s which is 
independent of the details of the initial condition (various vortex loops configurations were tried). Fig. [5] (top) shows 
a snapshot of such vortex tangle. 

Our vortex line densities agree with previous work; for example, taking T = 1.9 K we obtain 7 w 137 s/cm^ which 



compares well to 7 « 140 and 133 obtained in the numerical simulations of Adachi et al. [24| and in the experiments 
of Childers and Tough [3J| respectively. 

To analyze our results wc Fourier-transform the superfiuid velocity and compute the energy spectrum E'fc. If the 
turbulence is isotropic, Ek is defined by 

1 /■ 1 P°° 

E-yl-,<dV = l E,dk, (4) 

where V is volume, k = |k|, and k is the three-dimensional wavenumber. However, it is well known [2ll. |24{ that 
counterflow turbulence is flattened on the {y, z) plane perpendicular to the direction (x) of the heat flux. For example, 
if Lx, Ly and Lz are the vortex lengths (per unit volume) projected in the x, y and z-direction, at Vns ~ 1-25 cm/s 
we have L^/L = 0.34 < Ly/L = 0.55 = L^/L. It is therefore better to distinguish between parallel and perpendicular 
superfiuid energy spectra, E\\ and E^, calculated replacing v^ = vl^+vly+vl^ in Eq. (g]) with ?>v1^ and (3/2)(d^j^-|-'u^^) 
respectively. Fig. |3] (top) shows E±_ for various Vns plotted in the range kjj = 211/0 < k± < fcAj = 27r/A^ (where k± 
is the perpendicular wavevector); the vertical lines mark the wavenumbers ki = 2^/1 corresponding to the average 
intervortex spacing, i « L^^'^. It is apparent that the perpendicular energy spectrum E± has a broad peak in the 
mesoscales at intermediate wavenumbers kjj < k < kg. At larger k the spectrum follows the typical k~^ scaling 
of smooth isolated vortex lines as expected. The parallel spectrum E\\ vs k\\ (where fcy is the parallel wavevector) 
exhibits similar features, see Fig. |3] (bottom). Plotting E±^ and E\\ vs k rather than k± and k\\ yields similar results. 
We note that the counterflow energy spectrum, which we measure, is qualitatevcly similar to the spectrum shown by 
Nemirovskii, Tsubota and Araki [33 in their Fig. 2. 

Proceeding in analogy to what wc did for counterflow turbulence, we start from an arbitrary seeding initial condition, 
drive the vortex tangle with the synthetic turbulent flow of Eq. ([3]), and let L grow and saturate to a statistical steady 



state of turbulence which docs not depend on the initial condition (the time behaviour of L is similar to Fig. [T|). 
A snapshot of this mechanically-driven tangle is shown in Fig. [2] (bottom). We then compute the superfluid energy 
spectrum. In agreement with previous experimental [3, 0] and theoretical [3, [li| work, we find that the energy is 
concentrated at the largest scales, fc w fc^i, and that E^ ^ fc~^/^ for large k (see Fig.|3|). 

We conclude that there is a remarkable spectral difference between thermally-driven turbulence and mechanically- 
driven turbulence. Whereas in the former the turbulent kinetic energy is concentrated at intermediate length scales, 
in the latter most of the energy is at the largest scales, as in classical ordinary turbulence. 

An argument is often made in the literature that counterflow turbulence has only one characteristic length scale, 
the intervortex distance t. it is apparent from Fig. [3] that Ek does not have a sharp peak ai k k, ki = 2tt/(, (indicated 
by the vertical lines), but rather a broad maximum at smaller wavenumbers in the mesoscale region kjj < k < ki. 
The traditional argument, although quantitatively wrong, is thus qualitatively correct. 

IV. CURVATURE 

If we look carefully at the vortex tangles shown in Fig. [51 we notice that the thermally-driven tangle (top) contains 
relatively more closed loops, and the mechanically-driven tangle (bottom) contains relatively more long vortices which 
extend throughout the periodic computational domain. We sample the curvature C = |s"| along each vortex loop and 
construct the probability density function (PDF) of the mean curvature C of each distinct loop. Fig.[S]shows the result. 
We notice that mechanically-driven turbulence contains smaller curvatures (that is, larger radii of curvature R = l/C*) 
than thermally-driven turbulence; indeed, for the latter PDF(C) has a maximum at C « 250 in correspondence of 
the maximum of the energy spectrum shown in Fig. [3l 

As an additional numerical experiment, we compute the energy spectra of configurations of circular vortex rings 
placed randomly in the periodic box of size D as a function of the rings' radius R. We find that ii R > D (in which 
case rings are "folded" into broken arches by the periodic boundary conditions) most of the energy is concentrated 
at the largest length scales, whereas if _R < D the energy spectrum peaks at intermediate scales, in analogy with the 
counterflow spectrum. 

V. COHERENT STRUCTURES 

We also notice another difference between the two forms of turbulence. If we convolve the vortex filaments with 
a Gaussian kernel and define a smoothed vorticity field iVs (the details of the procedure are described in Ref. [101 ) > 
it becomes apparent - see Fig. |6]- that the thermally- induced tangle (sustained by the uniform v^^*) is essentially 
featureless, whereas the mechanically-induced tangle (sustained by the turbulent vjf'*) contains "vortical worms", or 
regions of concentrated vorticity - see Fig. [7l This result is consistent with the observation of "worms" in two other 
related turbulent flows: ordinary viscous turbulence |12| and pure superfluid turbulence at T = without the normal 
fluid |l3|; both flows satisfy the Kolmogorov k~^^^ scaling. 

It is known from previous work that if intense regions of normal fluid vorticity are imposed, such as Gaussian vortex 
tubes [3^ , ABC flows [331 o^ worms [S^, [S^ , these structures will induce (via the friction force) similar structures in 
the superfluid vortex lines. Our synthetic turbulent flow v^^*, although not completely featureless on its own, contains 
only weak vortex structures, much smaller [40[ than the vortical worms arising from direct numerical simulations of 
the Navier-Stokes equation. Therefore the observation of superfluid vortex bundles driven by the synthetic turbulent 
flow vfj^* of Eq. ([3|) must be an underestimate of the strength of these bundles. If we solved the Navier-Stokes equation 
for the normal fluid (rather than imposing v^^*), the normal fluid's worms would probably "imprint" vortex bundles 
in the superfluid, besides the bundles which arise naturally in the superfluid as a consequence of Euler dynamics ]1(^ . 

A tentative explanation of the observation that the vortex configuration is rather homogeneous for thermally-driven 
turbulence and inhomogeneous for mechanically-driven turbulence is that in the former (assuming, as we do, a uniform 
normal flow) the growth rate of the Donnelly-Glaberson (DG) instability (which transforms normal fluid's energy into 
superfluid vortex length) is the same everywhere, whereas in the latter it changes with time and space. 

The DG mechanism is the following [41|. If it is large enough, the component V of the normal fluid velocity along 
a vortex lines can destabilise a (helical) Kelvin wave of given wavenumber k. In this case, the Kelvin wave grows with 
amplitude A{t) = ^(O)e'^*, where A{0) is the initial amplitude of the helix and 

a{k) = a{kV ~ v'k^) (5) 

is the growth rate, v' = k,Ci/{A'k) « k, and £i = ln[l/(fcao)]. The growth of the Kelvin wave, however, may be 
interrupted by a vortex reconnection which "breaks" the vortex line. It is known that vortex reconnections play an 



essential role in the turbulence [21|, |42, |43[. Therefore, it is prudent to assess the effect of reconnections on the DG 
instability. 

Consider mechanically-driven turbulence in a statisticaly steady state at T = 1.9 K (a = 0.206) driven by the rms 
normal fluid velocity V « 0.93 cm/s, with average vortex length A « 11.5 cm, vortex line density L « 1.15 x 10'* cm^^, 
and intervortex spacing i « 9.3 x 10^'^ cm, and thermally-driven turbulence at the same temperature with V = 
0.75 cm/s, A « 11.88 cm, L w 1.19 x 10"' cm^^, and £ « 9.2 x 10^'^. The average number ^ of vortex reconnections 
per unit time is monitored during the numerical calculations; we obtain C « 4370 and 7386 s~* for mechanically and 
thermally driven turbulence, respectively (in reasonable agreement with the estimate ( ~ (2/3) kL^/^ w 9500 s~*, for 
a homogeneous isotropic tangle, of Barenghi & Samuels |44|V 

The mode which undergoes the most rapid DG instability has wavenumber fcmax = ^/(2i^') and growth rate 
CTmax — ctV'^ /{Av'), Corresponding to the length scale dmax — 27r/fcmax. In both mechanically-driven and thermally- 
driven cases this length scale (dmax = 0.015 and 0.017 cm respectively) is larger than than the average distance 
between vortices £, and so not suitable for our analysis. We therefore perform an analysis for Kelvin waves with a 
wavelength and a wavenumber equal to i and kg = 27r/£, respectively; we assume that such waves are the lowest 
frequency waves in our system. The growth rate of such a waves is (Tdq = cr^ki), where <7{k) is defined in Eq. ([5|). 

The reconnection rate C computed during the simulations is a statistical property of the vortex tangle as a whole. 
However we can compute a reconnection frequency for a wavenumber and wave amplitude by scaling the total re- 
connection rate by the fraction of the total vortex length that a given wavelength takes. In such a manner we 
define 



ar{A) = ^ f [1 + (Aki)^ cos'^ikix)] '^' dx. 
^ Jo 



(6) 



Fig.[5]shows the ratio ooGJOr vs wave amplitude A for the two simulations described above. We note a large contrast 
in the behaviour of the ratio of these two timescales when comparing the mechanically and thermally driven cases. 
For the latter one would estimate that the amplitude of perturbations along the vortices can grow to approximately 
the intervortex spacing before reconnections dominate the behaviour of the tangle. However, in the mechanically 
driven case odg ~ ^r for A ~ D/3 so that the large amplitude perturbations are able to grow, before reconnections 
randomise the tangle and introduce topological changes. Therefore the difference in the balance between these two 
competing timescales is likely to be partially responsible for the differences in the nature of the two turbulent systems. 

VI. VORTEX RECONNECTIONS 



The existence of superfluid vortex bundles [36|-|39| , their dynamics [45| and their particular significance at very low 
temperatures |46l. |47|| have been discussed in the literature, but so far there is no clear experimental evidence for them. 
It has been argued that the presence of bundles of locally almost parallel vortices (which we have demonstrated in 
the previous section for mechanically- induced turbulence) leads to a suppression of vortex reconnections J46l |47{ . 

In the vortex filament model, vortex reconnections are performed algorithmically; the details are described in 
Ref. |3l[ . Within the approximation, intrinsic to the model, it is instructive to study the distribution of the angles 9 
between reconnecting vortex lines at the level of discretization which we use (which is necessarily much larger than cq). 
The normalised distribution of values of 9, PDF(6'), is shown in Fig.[Hl the solid black line with black circles refers to 
mechanically-driven turbulence, and the solid red line with red squares to thermally driven turbulence. It is apparent 
that in thermally-induced turbulence the majority of vortex reconnections take place between vortex filaments which 
are nearly anti- parallel (0 w tt), whereas in mechanically-driven turbulence most reconnections are between vortices 
which are nearly parallel (9 < tt/2). Our results confirm that indeed the presence of organised bundles of vortices 
changes the typical geometry of reconnections. 

To check the temperature dependence of the results we repeat our calculations at higher temperature, T ~ 2.1 K. 
At this temperature the friction coefficients are larger {a = 1.21 and a' = —0.3883), therefore a more intense vortex 
tangle is generated at the same value of the drive; moreover, short Kelvin waves are damped out more quickly. We 
check that at the higher temperature T = 2.1 i^ the same differences between thermally-driven and mechanically- 
driven turbulence are present, which we have described in the previous sections for T = 1.9 K in terms of energy 
spectrum, curvature and coherence structures. Fig.[3]shows that, qualitatively, the distribution of reconnecting angles 
is also temperature independent (the black solid line, which refers to mechanically-driven turbulence peaks at small 
9, the solid red line which refers to thermally-driven turbulence peaks at large 9). 

This result could be exploited to look for experimental evidence of superfluid vortex bundles in the following way. 
Using solid hydrogen tracer particles to visualise the vortex lines, Paoletti, Fisher and Lathrop et al. [4^ determined 
that the the minimum distance 6(t) between vortex lines before and after a reconnection scales as 



5{t) ^ A{K\t - to\y^^{l + c\t ^ to\), (7) 

where to is the time at which the reconnection takes place, with fitting coefficients A ?» 1.2 and c w 0. We proceed 
in this way, monitoring vortex reconnections in our numerical calculations. Fig. [TU] shows the probability density 
functions of our fitting parameters A and c obtained for 1107 reconnections in thermally-driven turbulence (average 
values (A) ~ 2.6 and (c) = 1.6) and 879 reconnections in mechanically-driven turbulence (average values (A) ~ 1.8 
and (c) — 0.7 s""'^). Our fitting coefficients thus agree fairly well with the experimental findings of Ref. [43 and with 
the numerical results of Tsubota and Adachi [49[ (A w 3 and c « s~^). 

Fig. [To] shows that the distribution of values of A is different for thermally-driven and mechanically-driven turbu- 
lence. The effect must arise from the different distributions of curvature and reconnecting angles 9 for vortex bundles, 
which are present only in mechanically-driven turbulence. This is confirmed by Fig. 111! which displays scatter plots 
of the fitting parameters A (top) and c (middle). Fig. [TTjalso shows the angular dependence of the mean curvature 
Cr (bottom) of the reconnecting vortex segments. 

It is clear that the curvature of the filaments after a reconnection is dependent on the angle of the reconnection. From 
inspection of the local induction approximation [23| we would expect larger velocities (and thus A) with increased 
curvature, as we observe in the numerical simulations. These result suggests a possible experimental strategy to 
establish the existence of vortex bundles based on the careful analysis of the reconnection fitting parameter A. 

VII. CONCLUSIONS 

In conclusion, we have addressed for the first time the question of the energy spectrum of thermally-induced 
counterflow turbulence, and found that it is unlike the spectrum of turbulence generated mechanically. More in 
general, we have found that the two forms of quantum turbulence which can be generated in superfluid helium are 
quite different. Counterflow turbulence, driven thermally by a constant normal fluid velocity, is uniform in physical 
space and the energy spectrum is concentrated at intermediate wavenumbers k. On the contrary, quantum turbulence 
driven mechanically by a turbulent normal fluid contains regions of concentrated coherent vorticity and vortex lines 
with larger radii of curvature; the energy is concentrated at the largest scales, exhibiting the same k~^'^ scaling of 
ordinary turbulence which suggests the presence of an energy cascade. Our results prove that counterflow turbulence, 
a form of disordered heat transfer unique to liquid helium, lacks the multitude of interacting length scales which is 
perhaps the main property of ordinary turbulence. Vortex reconnections are affected by the presence of bundles of 
almost parallel vortices, suggesting an experimental technique to detect these bundles based on monitoring the vortex 
separation after reconnections. 
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FIG. 1. (Color online) Thermally-induced turbulence. The evolution of the vortex line density L (cm ^) vs time t (s) at 
counterflow velocities (from top to bottom) Vna ~ 1-25 cm/s (black), 1.0 cm/s (red), 0.75 cm/s (blue), and 0.55 cm/s (green). 





FIG. 2. Snapshots of vortex tangles (y,z projections). Top: thermally-driven by counterflow (u„ 
12000 cm"^); bottom: mechanically-driven {Re = 208, L « 12000 cm"^) 
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FIG. 3. (Color online) Counterflow turbulence. Top: the perpendicular energy spectrum E_i_{k±) (arbitrary units) vs 
wavenumber kj_ (cm^^). Bottom: the parallel energy spectrum i5||(fc||) (arbitrary units) vs wavenumber fcy (cm^^). The 
vertical lines mark ke at increasing v„s from right to left. 
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FIG. 4. (Color online) Mechanically-induced turbulence. Energy spectrum Ek (arbitrary units) vs wavenumber k (cm ^) 
of vortex tangle driven by the synthetic turbulent flow of Eq. ((3)1 with M = 188 modes. The vertical dashed blue line 
marks ke. The dashed red line shows the fc~^'^ Kolmogorov scaling. The effective Reynolds number of the normal fluid is 
Re = (fcAz/fci)*'"^ — 208, deflned by the condition that the dissipation time equals the eddy turnover time at k^j. 
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FIG. 5. Probability density function (PDF) of the mean curvature per vortex loop C. Solid black line: mechanically-driven 
turbulence; dashed red line: thermally-driven turbulence. Notice the larger curvatures present in thermally-driven turbulence. 
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FIG. 6. (Color online) Smoothed vorticity u! sustained by a constant vj^^* (thermally-driven turbulence). 




FIG. 7. (Color online) Smoothed vorticity oj sustained by a turbulent v^^ 
intense vortical regions compared to Fig. |6] which is plotted on the same scale. 



(mechanically-driven turbulence). Notice the 
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FIG. 8. (Color online) Plot of aDcl^r (ratio of Donnelly-Glaberson and vortex reconnection frequencies) vs wave amplitude 
A (cm), for thermally (dashed, red line) and mechanically (solid line) driven turbulence; the (blue) dot-dashed line represents 
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FIG. 9. (Color online) The probability density function (PDF) of the angle between reconnecting vortices, 6. Thermally-driven 
turbulence: at T = 2.1 K (dashed red line) and at T = 1.9 K (solid red squares); mechanically-driven turbulence: at T = 2.1 K 
(solid black line) and at T = 1.9 K (solid black line with solid black circles). Note that for thermally-driven turbulence the 
distribution peaks at large 0, whereas for mechanically-driven turbulence it peaks at small 0. 
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FIG. 10. (color online) Probability density functions (PDF) of the fitting parameters A (top) and c (s ^) (bottom) of Eq. [7] 
Solid black line: mechanically-driven turbulence; dashed red line: thermally-driven turbulence. 
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FIG. 11. (color online) Scatter plots of the fitting parameters A (top) and c (middle) of Eq. [7]vs the angle, 0, between the 
reconnecting vortices. The bottom figure shows the mean curvature Cr of the reconnecting vortex segments vs the angle 9. 
Solid black points: mechanically-driven turbulence; open red circles: thermally-driven turbulence. 



